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Abstract 

The linear and nonlinear pulse propagations in lifetime-broadened three-state media with spon- 
taneously generated coherence (SGC) are investigated theoretically. Three generic systems of V-, 
A-, and H-type level configurations are considered and compared. It is shown that in linear prop- 
agation regime the SGC in the V-type system can result in a significant change of dispersion and 
absorption and may be used to completely eliminate the absorption and largely reduce the group 
velocity of a probe field. However, the SGC has no effect on the dispersion and absorption of the 
A- and H-type systems. In nonlinear propagation regime, the SGC displays different infiuences on 
Kerr nonlinearity for different systems. Specifically, it can enhance the Kerr nonlinearity of the 
V-type system whereas weaken the Kerr nonlinearity of the A- type system. Using the SGC, stable 
optical solitons with ultraslow propagating velocity and very low pump power can be produced in 
the V-type system by exploiting only one laser field in the system. 



PACS numbers: 42.65.Tg, 05.45.Yv 



I. INTRODUCTION 



The propagation of linear and nonlinear light pulses in coherent atomic media has been 



an irajjortant subject of many recent studies, such as formation o: 
tons [1] , group- velocity reduction ^ , storage and retrieval of light 



simultons and adiaba- 
y] , and single-photon 



pulse propagation jsj. Among various coherent preparation techniques studied so far, elec- 
tromagnetic induced transparency (EIT) is perhaps the most extensively investigated one 
because of its diverse practical applications {g]. Due to the quantum interference induced 
by a coupling laser field, a probe laser field propagating in an EIT medium can avoid a 
large absorption even when it is tuned to a strong one-photon resonance. EIT media also 
exhibit drastic change of dispersion and giant enhancement of Kerr nonlinearity, resulting 
in various nonlinear optical phenomena such as high-efficient frequency conversions [3, Isl, 
temporal and spatial optical solitons at very low light level js iSl. As a powerful technique, 
EIT has been used to study the property of Rydberg atomic ensembles very recently [l3|. 

In an EIT medium, it is crucial to have at least two laser fields as they are necessary to be 
used to create atomic coherence. Besides the EIT technique, an atomic coherence can also 
be created by the quantum interference between two spontaneous emission channels without 
using any coupling laser field, which is called the spontaneously generated coherence (SGC). 
In recent years, much attention has been paid to the study on SGC and related topics, includ- 

16|, coherent population trapping [l7|, spectral narrowing 
fiuorescence squeezing [22], giant self-phase modulation 



ing lasing without inversion [14 
and fiuorescence quenching [18 



23|. ground-state quantum beats '24|, cavity- mode entanglement 25], electromagnetically 



induced grating 



26|, and so on. 



Although a large amount of research activities have been made on SGC, most of them are, 
however, focused on the static property of various systems, and only a few works dedicate 
to the study of pulse propagation in SGC media. Here we mention the work by Paspalakis 
et al. [27^ who studied the pulse propagation in a four-level system, where a ground state 
is coupled to two closely spaced excited states by a laser field with both excited states 
decaying into a common continuum. Strong, short laser pulses with adiabaton-like property 
were observed by using numerical simulations. 

In this article, we investigate, both analytical and numerically, the linear and nonlinear 
pulse propagations in lifetime-broadened three-state media with SGC. In stead of adiabatons. 
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we consider breather-like nonlinear excitations without using any adiabatic approximation. 
Our work includes two aspects: (i) We consider three generic systems of V-, A-, and H-type 
level configurations, and compare their linear propagating property with the SGC effect 
being taken into account. We show that the SGC in the V-type system can result in a 
significant change of dispersion and absorption and may be used to completely eliminate the 
absorption and largely reduce the group velocity of probe field. However, the SGC has no 
effect on the dispersion and absorption of the A- and H-type systems, (ii) We demonstrate 
that in nonlinear propagation regime, the SGC has different influences on Kerr nonlinearity 
for different systems. In particular, it can enhance the Kerr nonlinearity of the V-type system 
but weaken the Kerr nonlinearity of the A-type system. By using the SGC, stable optical 
solitons with ultraslow propagating velocity and very low pump power can be produced 



br generating ultraslow optical solitons 
9l-ll2| because the suppression of optical 



in the V-type system. We stress that the scheme 
presented here is very different from those in Refs. 
absorption and the reduction of group velocity are not contributed by an additional control 
fleld but by the SGC and hence only a single laser fleld is needed. 

The work reported here is arranged as follows. In Sec. II, three-level models of V-, A-, 
and S-type conflgurations with SGC are introduced. Linear pulse propagations are discussed 
and dispersion and absorption properties for three different systems are analyzed. In Sec. 
Ill, the Kerr nonlinearity of the V- and A-type systems are analyzed, and ultraslow optical 
solitons at very low light level are obtained. In Sec. IV, a discussion of open system is made. 
The last section contains a summary of our main results. 



II. MODELS AND PULSE PROPAGATION IN LINEAR REGIME 

For comparison and also for completeness, we investigate three generic three-state systems 
of V-, A-, and H-type level conflgurations, which are considered separately in the following. 



A. V-type system 

We flrst consider a three-level V-type atomic system, as shown in Fig. [Il^a), in which 
two closely spaced excited states |2) and |3) decay simultaneously into the ground state 
|1) by the spontaneous emission with decay rates r2 and T^, respectively. The quantum 
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FIG. 1: (color online) Energy level diagrams and excitation schemes of lifetime-broadened three- 
level systems with SGC. (a): V-type system; (b): A-type system; (c): H-type system. \j) (j = 
1, 2, 3) are atomic bare states, Ep is weak probe laser field, A and 5 are detunings, and Tj (j = 1, 2) 
are decay rates of relevant states, (d): The definition of the alignment angles (^i, 62) of the dipole 
matrix elements (di2, dis) related to the unit polarization vector of the probe field. 9 = 9i + 62 
is the angle between di2 and dis. 

interference between the two decay channels (from |2) to |1) and |3) to |1)) resuhs in the 
SGC of the system [28|. A weak, pulsed probe field (with duration tq) of center frequency 



oOp and wavevector kp, i.e., 



Ep(r,t) = ep£p{r,t)e 



i(kp-r— Wpi) 



+ C.C, 



(1) 



couples the ground state |1) to the excited states |2) and |3), where and £p{r,t) are the 
unit polarization vector and envelope function of the probe field, respectively. 

Under electric-dipole, rotating-wave, and Weisskopf-Wigner approximations, the equa- 



tions of motion for the density matrix governing atomic dynamics are 



P22 
P33 
P21 

P31 
P32 



-r2P22 + i^pPu - i^*pP21 - V- 



-r3P33 + ip^pPis - ip^pPsi 

£2" 

2 

£3 
2 

r2 + r3 



V- 



2 



(P23 + P32), 
(P23 + P32), 



i2A 



P21 + i^piPii - P22) - ip^pP23 - V—K-^p3i, 



p3l + ipQpipn - P33) - ^fipP32 - V ^^^^^ 21, 



^ P32 - i%P31 + ip^pPl2 - V—^Y^{P22 + P33 



(2a) 
(2b) 
(2c) 

(2d) 

(2e) 



with pii + P22 + P33 = 1- Here Qp = ep ■ di2£p/h is half Rabi frequency of the probe field 
with dij = (^|d|j) being the density-matrix elements related to states \i) and A = 
{E^ — E2) I {2h) is half frequency difference between |2) and |3), and 6 = Up— {E3 + E2) / {2h) 
is one-photon detuning [see Fig. IH^a)]. The cross coupling term contributed by the SGC 
effect is manifested by the factor rj^/^2^/2, with rj = di2 ■ di3/|di2| |di3| = cos6' denoting 
the alignment of two dipole matrix elements di2 and di3, where 6 is the misalignment angle 
between di2 and di3. If di2 and di3 are parallel (i.e. 6 = 0), one has rj = 1, the system 
exhibits maximum SGC; if di2 and di3 are perpendicular (i.e. 9 = 7r/2), one has rj = 0, the 



system displays no SGC. p = |ep ■ d 



131 



|e„ ■ d 



121 



|di3| cos^^i/[|di2| cos6'2], where 9i (6*2) is 



the misalignment angle between di2 (di3) and Bp. In the following, we assume |d 
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and a particular case can be found that 6 is equally partitioned by e^, i.e. 61 ^ 62 



|di2| 

e/2, 

as did by Wan et al. [26|. In such case, we have p ~ 1 and r] is still kept a free parameter 
with its value taking between —1 and 1 [see Fig. [11(d)]. 

The equation of motion for the probe-field Rabi frequency Qp can be obtained by using 
Maxwell equation. Under slowly-varying envelope approximation, it reads 



d ^Id 
dz c dt 



^p + k{p2i + pp3i) = 0, (3) 

where k = AfaUJp\ep ■ di2p/(2eoc/l) with Afa being the atomic concentration. For simplicity, 
we assume in the following that r2 ~ r3 = F. 

The linear optical response of the system can be obtained by solving the Maxwell-Bloch 
(MB) Eqs. ([2]) and ([3]). Assuming Qp is a small quantity, pn ~ 1, and P21, P31, and Qp are 
proportional to exp[i{Kz — ut)], we obtain the linear dispersion relation 

CO + ds — ipriT/2 ^ p^iu + (^2) — ipriT/2 



D u 



D{u) 



(4) 
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FIG. 2: (color online) (a), (b): Re(-ftr) and Im(if) as functions of w/F with maximum SGC (i.e. 
r/ = 1; the solid line) and without SGC (i.e. r/ = 0; the dashed line), (c), (d): Im(ii'o)|5=o and 
(V"g/c)|5=o as functions of ??. 

where D(u;) = (w + rf2)(a; + 4) + rfV^jl with cis = A + 5 + iir/2, and 4 = -A + 5 + ir/2. 
By Taylor expanding K(uj) around = [29], we obtain K(uj) = + -ft'iw + \K2iJ^ + • • • , 



with the expansion coefficients Kj 
p'^d2 + ds- ipr]T 



Kq = -K 
K2 



4^3 + r/'rV4 

1 1 



+ K 



[9^-ir(a;)/9u;^-]Uo (j = 0, 1, 2,- 



(^^(^2 + C?3 - ^pr^r) {d2 + (^3 



exphcitly given by 
(5a) 
(5b) 



(l+p2)(rf2 + c/3) p^d2 + ds 

2k—-— „„„ , + 2k- 



ipr]T {p'^d2 + dj,- ipr]T) (^2 + d^ 



2k 



-(5c) 



(^24 + r/'rV4)2 {d2ds + r]^TyAy {d2ds + v^T^ 4)^ 

Here, ii'o = Re(-K'o) + ilm{KQ) gives the phase shift per unit length and absorption coeffi- 
cient, Ki determines the group velocity Vg (= 1/Ki), and K2 represents the group- velocity 
dispersion. 

Shown in Shown in Fig. [2]^a) and Fig. EJ^b) are respectively the real part and imaginary 
part of K as functions of u, which characterize the dispersion and absorption of the system. 
The solid lines in the figure are for the case with the maximum SGC (i.e. rj = 1), whereas 
the dashed lines are for the case without the SGC (i.e. 77 = 0). System parameters are 
chosen as k = 1.0 x 10^ cm-is~\ F = 1.0 x 10^ s-\ A = 1.0 x 10^ 5 = 0, and p = 1. We 
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see that in the region around u = the probe-field displays a drastic change of dispersion 
(and hence a drastic reduction of group velocity) (panel (a) ) and a large suppression of 
absorption (panel (b)). Obviously, the reduction of group- velocity and the suppression of 
absorption with the SGC are much more significant than those without the SGC These can 
be seen also by the expression of the group velocity and the absorption coefficient at 5 = 0: 

ri A2-(i-n)2rV4 

Mi^o)i.^o-^,^[^^i;^, (6b) 

which are respectively shown in the panels (c) and (d) of Fig. |2l Notice that the group 
velocity of the probe field can be lowered by increasing the SGC effect (panel (c) ). li rj = 1 
one has Im(i^o)|5=o = 0, i.e. the absorption of the probe field can be completely eliminated 
by means of the SGC. But if the SGC is absent (i.e. t] = 0), lm.{KQ)\s=o is not only non-zero 
but also with a large positive value, and hence the probe-field absorption is quite significant 
(panel (d) ). Consequently, in the absence of the SGC the probe filed can not propagate to 
a long distance. 

It is instructive to discuss the SGC effect on the probe-field absorption in more details. 
From Fig. [2](b), we see that the probe-field absorption is suppressed around u = in both 
cases with and without the SGC, resulting in the appearance of transparency windows in 
the absorption spectrum. However, the depth and width of the transparency windows are 
quite different. The transparency window with the SGC [r] = 1; the solid line) is much more 
deeper and wider than that without the SGC (r/ = 0; the dashed line). 

The difference for the absorption spectra for the above two cases can be understood more 
clearly as follows. For illustration, we split the imaginary part of K into the form 



Im{K) = K 
+ 



F / 1 1 

+ 



_2 V(w + 5-t;)2 + rV4 (u; + 5 + t;)2 + rV4 
rjT f u + 6 — V UJ + 6 + V 



(7) 



with V = a/A^ — rj'^T'^/A. The first two terms on the right hand side (RHS) of Eq. (JTj) 
correspond to two resonances of the excited states |2) and |3), shown by the dashed line of 
Fig. [3l The dip with a non-zero minimum in the dashed line can be interpreted as a gap 
between the two resonances, which is a typical character of Autler-Townes (AT) splitting 
3o| . The next two terms (interference terms) provide a gain (absorption) when rj > 




FIG. 3: (color online) The spectra of AT-splitting (dashed line), SGC (dotted line), and total 
absorption (solid line) as functions of lo/T with 5 = 0. 

[t] < 0), shown in Fig. |3] by the dotted-dashed hne. Such terms come from the quantum 
interference effect between two spontaneous emission channels (i.e. from |2) to |1) and from 
1 3) to |1)), a typical character of SGC. The sum of the four terms on the RHS of Eq. ([7]) is 
the total absorption of the probe field, shown by the solid line of Fig. [31 From these results 
we clearly see that it is the joint contribution from the AT splitting and the SGC effect 
that make the absorption of the probe field vanish, and hence providing the possibility for 
a long-distance transmission of the probe field in both the linear and nonlinear propagation 
regimes. 



B. A-type system 

We now consider a three-level A-type system, in which a weak probe field with the form 
([I]) couples the excited state |3) to the two ground states |1) and |2), with corresponding 
spontaneous decay rates Fi and F2, respectively (see Fig. [11(b)). A SGC occurs by the 
quantum interference between two spontaneous decay channels from |3) to |1) and from |3) 



to |2) |3l|,l32|. 
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The atomic dynamics of the system is described by the density matrix equations 



Pll 
P22 
P21 

P31 
P32 



riP33 + ip%Pzi - ip^ppi3, 

r2P33 + ^^pP32 - ^^pP23, 



-i2Ap2i + i^*pp3i - ip^pp23 + r]V^1^2p33, 



-% (A + b) 
i(A-5) - 



2 

ri + r2^ 



P31 + ip^lp{pn - P33) + «J^pP2i, 

P32 + i^p{p22 - P33) + ip^Pu, 



(8a) 
(8b) 
(8c) 
(8d) 

(8e) 



with + P22 + P33 = 1. Here Qp = Bp ■ di^Sp/h is half Rabi frequency of the probe field, 
A = {E2 — Ei)/{2h) is half frequency difference of two closely spaced ground-state levels, 
and 5 = — Ei — E2)/{2h) — Up is one-photon detuning (see Fig. [T]^b) ). The SGC effect 
is still described by the factor rjy/ViV2/'2,, with rj = dis ■ d23/|di3| |d23|. 



The equation of motion for Qp is 



+ ^. + '^(^/^31+P32)=0, 



(9) 



where p = \ep ■ di3|/|ep ■ d23| = cos 611/ cos 612 (|di3| ~ |d23|) and k = AfaUp\ep ■ d23p/(2eoc^). 
For simplicity, we assume in the following that Fi ^ r2 = F. 

From the MB Eqs. ([8]) and Q it is easy to obtain the linear dispersion relation of the 
system, which reads 



c 



K 



P 



2 Pll 



+ 



(0) 
P22 



where di = A + 6 — iT, ^2 



u — di u — d2 
A + 5 - «F, pS? = Mi|V(Mil' 



(10) 



(0) 



^^1^2! ), and p^2 



P^M2p/(MiP + P^\d2?)- As before, we Taylor expand Kioj) at w = 0, i.e. Kioj) = Kq + 
Kiu + i/sTsw^ + ■ ■ ■ . The expansion coefficients Kj = [d^ K{u)/duj^]\uj=o (j = 0, 1, 2,---) 



are given by Kq = n{p'^pf}/di + P22 /d2), Ki = 1/c + K{p'^Pu/dl + P22 /^l), and K2 = 

2«:(pVSiVrf? + pgV4)- 

Equation ( ITOj) consists of two simple Lorentzian terms, corresponding to the resonances 

between |3) and |1) and between |3) and |2), respectively. Thus only an AT splitting in 

the absorption spectrum of the probe field occurs. This tells us in the A-type system the 

SGC has no effect on the linear dispersion and absorption spectra of the system. Shown in 

Shown in Fig. Hl^a) and Fig. Ht^b) are Re(i^) and Im(i^) as functions of u. We see that a 

dip appears in Im(i^) (panel (b) ) and, at the same time, the dispersion changes rapidly in 



.(0) 



2.(0) 



.(0) 




FIG. 4: (color online) Jie{K) (panel (a) ) and Im(K) (panel (b) ) as functions of u/T with A = F 
(solid line) and A = 2F (dashed line). 

the region around u = ((panel (a) ). In addition, Re(i^) and Im(i^') can be adjusted by 
changing the value of A. The solid (dashed) line in the figure is for A = F (A = 2F). We 
see that although lm{K) can be lowered by increasing A, it has always a large value. Thus 
in the present A-type system one can not acquire a small absorption for the probe field. We 
must also point out that the value of A can not be adjusted to be too large. This is because 
for a large A the SGC effect disappears. 



C. H-type system 



We finally consider a three- level S-type system, as shown in Fig. [I](c), in which a weak 



2)_ and simultaneously 
States 1 2) and |3) 



probe field couples the ground state |1) to the intermediate state 
couples the state |2) to the excited state |3), as suggested in Ref. 
decay to |1) and |2) with decay rates Fi and F2, respectively. The MB equations governing 
the evolution of the atoms and the electric field are given by 



P33 = -rsPsS + ip^pp23 - W%P32 

P21 = ^-^(^ - P21 + i^pipn - P22) + ip^pPsi + vV"^^ 

P31 + ip^pP21 - i^pP32 

r2 + ri 



3P32, 



Psi 

P32 



.,2(A + 5)-^ 



-i(2A + 5)- 



P32 + ipVtp{p22 - P33) - i%p3l, 



M ^ + ) + '«(P2i + PP32) = 0, 



(lla) 
(lib) 
(11c) 

(lid) 

(lie) 

(llf) 
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with pii + P22 + P33 = 1, where Qp = \ep ■ di2\Sp/h is half Rabi frequency of the probe 
field, A = {E^ — 2E2 + Ei)/{2h) is half frequency difference of the transitions |1) o |2) and 
|2) ^ |3), and 6 = {E2 — Ei)/h — cOp is one-photon detuning [see Fig. HJ^c)]. The last term 
on the RHS of Eq. flllcp comes from the SGC effect, with rj = di2 ■ d23/|di2| |d23| = cos6'. 
In Eq. ([n]), p=\ep- daal/le^ ■ di2|. 

Using the MB Eq. (11 II) we obtain the linear dispersion relation of the system 

K{u) = -- K / ^ , ■ (12) 

From this formula, we see that: (i) The SGC in the S-type system also does not change 
the linear dispersion and absorption spectra of the probe field; (ii) No dip structure in the 
absorption spectrum is observed because the system does not possess any AT splitting. 
Because the peak of the absorption spectrum locates at the center frequency of the probe 
pulse, a long-distance wave propagation of the probe pulse is not possible in the system. 
This is not interesting for wave-propagation problem and hence we discard this model in the 
following discussion. 

III. PULSE PROPAGATION IN NONLINEAR REGIME 
A. V-type system 

Kerr nonlinearity is essential for most nonlinear optical processes. It can be largely 
enhanced in resonant optical media, but usually a serious optical absorption is also accom- 
panied simultaneously. However, here we show that by the joint action of the AT splitting 
and the SGC effect, the Kerr nonlinearity in the V-type system can be enhanced greatly 
with the optical absorption eliminated greatly. 

The probe-field susceptibility for the V-type system is defined as 

_ A/;|ep-di2pP21+PP31 _ (1) (3) IP ,2 /.gx 
Ap — ^ — Ap ^ Ap l^pl 1 \^^) 

where Xp^^ ^ind Xp^^ are linear and third-order susceptibilities, respectively. The real part of 

(3) (3) 

Xp contributes to the Kerr nonlinearity while the imaginary part of Xp contributes to the 
nonlinear absorption or gain of the system. The explicit expressions of Xp'^ and Xp^ can be 
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FIG. 5: (color online) (a) and (b): Re(x*^^)) and lm.{x^^^) as functions of rj for the V-type system, 
respectively, (c) and (d): Re(x'-^^) and Im(x'^^'') as functions of for the A-type system, respectively. 



obtained by solving Eq. ([2]) under steady-state approximation, which reads 

A/'alep ■ di2pp^(i2 + 4 - ipV^ 



y(i) 

A.p 



Afa\ep- di2|^ d2A + 0^35 - iriTC/2 



(14a) 
(14b) 



60^3 c/arfs + ??2rV4 ' 

with the explicit expressions of A, B, and C being given in Appendix A. Shown in Fig. |5]^a) 
and Fig. |5]^b) are Ile(x^^^) and lm{x^^^) as functions of r], respectively. System parameters 
are taken as Afa = 10^ cm'^, dis ^ ^23 = 2.5 x lO'^^ cmC, F = 1.0 x 10^ s-\ A = 5.5 x 10^ 
s^^, 6 = 1.0 X 10'' s"^, and p = 1. From the figure we see that |Re(x^^^)| grows as rj mcreases, 
whereas |Im(x'^'^^)| is reduced as r] increases. The condition |Re(x^^'')| ^ holds in 

the whole range of rj. Thus, the SGC effect enhances the Kerr nonlinearity of the system 



significantly. In addition, Re(x^^'') has an order of IC 



cm 



-1 o2 



s , i.e. it is 10 times larger 



than that of conventional nonlinear optical media 34| . 

The Kerr nonlinearity enhancement obtained above can be used to balance the dispersion 
of the system and hence to obtain a lossless and distortionless optical pulse propagation in 
nonlinear regime. To this end we apply the standard multiple-scale method |9|, which is 
beyond the steady-state and adiabatic approximations, to solve Eqs. ([2]) and (|3]). We make 
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the following asymptotic expansion pjj = 6ij + Zl^i ^"Pji"* = 1^ 2, 3), pij = 

= 1, 2, 3; i 7^ j), and fip = X^^i where e is a small parameter characterizing the 

small population depletion of the ground state. To obtain a divergence-free expansion, all 
quantities on the right hand side of the expansion are considered as functions of the multi- 
scale variables zi = e^z {I = 0, 1, 2) and ti = eH {I = 0, 1). Substituting the expansion and the 
multi-scale variables into Eqs. ([2]) and ([3]), we obtain a chain of linear, but inhomogeneous 
equations which can be solved order by order. 

At the leading order, we get the linear solution Vt^p'^ = F exp{i[K{uj)zo — ooto]} and the 
dispersion relation, given by Eq. (jlj). At the second order, a divergence-free condition 
requires dF/dzi + {l/Vg)dF/dti = 0. Here F is a yet to be determined envelope function 
depending on the slow variables ti, Zi and Z2. At the third order, we obtain the nonlinear 
equation for F: 

^sr^-^^'WeM~-p.,)F\F\^^,, (15) 

where /3 = e~^/3 with /3 = 21m(ii'o) and 



d^A + d^B - i7^TC/2 

= —K, — — TT- . (16 



After returning to original variables, Eq. (ITSl) becomes 



where t = t — z/Vg and U = eFe'^^^^"^. Equation ( !T7|) usually has complex coefficients due 

to the resonant character of the system. However, as we shall show below, under the joint 

action of the SGC and the AT splitting, practical set of system parameters can be found to 

make the imaginary part of the coefficients be much smaller than their real part. Then Eq. 

f lT7|l can be approximated as a nonlinear Schrodinger (NLS) equation, which allows soliton 

solutions being able to propagate for a rather long distance without significant attenuation 

and distortion. The dimensionless form of Eq. (|T7]) is 

du d'^u , , . 

?Tr + + 2m M = «2/iM, 18) 
OS oa^ 

where s = —z/{2Lci), cr = t/tq, and u = U/Uq, jj, = Ld/La, with Ljj = Tq/\K2\ being 
the characteristic dispersion length. La = the characteristic absorption length, and 



Uq = {l/To)y \K2/W\ the characteristic Rabi frequency of the probe field. The symbol tilde 
denotes the real part of the corresponding quantity. 
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If Ld is much less than La (i.e., /x <^ 1, which is the case in the presence of the SGC), 
the term on the right-hand side of Eq. f|T8l) can be treated as a small perturbation and can 
be neglected at the first order. Hence Eq. (fTSl) reduces to the standard NLS equation, which 
is completely integrable and allows multi-soliton solutions. After returning to the original 
variables, a single soliton solution of the NLS equation corresponds to 





exp 



10Z — l- 



2L 



D 



(19) 



We now present a practical numerical example to support the above results. Consider a 
cold alkali atomic gas, for which the system parameters can be taken as ~ ^34 = 2.5 x 
10-27 cm C, K = 1.0 X 10^ cm-i s'S Eg ^ Eg = 1.0 x 10^ s'S A = 5.5 x 10^ s-\ 5 = 1.0 x 10^ 
s-\ r/ = 1, and p = 1. Then we obtain = 6.83 + ^0.23 cm-\ Ki = (0.73 - i0.05) x 10^^ 
cm-^s, K2 = (0.14 + Z0.06) x lO^^^ cm^^s^, and W = (1.00 + z0.03) x 10^^^ cm^^s^. Notice 
that the imaginary parts of these quantities are indeed much smaller than their real parts. 
The reason is that the SGC plays an important role to make the optical absorption of the 
system be largely eliminated. Taking tq = 0.8 x 10"'' s, we have = 0.45 cm, Uq = 1.5 x 10^ 
s~^, and La_ = 4.3 cm (// ^ 0.1). Then the propagating velocity of the soliton is given by 



Vg = 4.6 X 10"^ c. 



(20) 



i.e., the soliton formed in the V-type system travels indeed with an ultraslow velocity. We 
must stress that the present scheme for generating the ultraslow optical soliton needs only 
one laser field, which is different from the EIT scheme, where at least two laser fields are 



required 



9-12 



With the above parameters, it is easy to estimate the peak power of the ultraslow optical 
soliton by using the Foyntings vector, which reads 



Ppeak = 0.65/iW^, 



(21) 



where the cross-section area of the probe laser beam is taken to be vr x 10"'^ cm^. Thus, we 
see that very low input power is needed for generating the ultraslow optical soliton due to 
the enhancement of Kerr nonlinearity. 

We have also studied the stability of the ultraslow optical soliton presented above by 
using numerical simulations. Shown in Eig. [6]^a) is the wave shape of {Qp/Uol"^ as a function 
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(a) (b) 
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FIG. 6: (color online) (a): Waveshape of |Op/C/oP as a function of z/L£) and t/tq. The solution 
is numerically obtained from Eq. (I17p with full complex coefficients taken into account, (b): 
Collision between two ultraslow solitons. (c): Waveshape of |r2p/C/oP as a function of t/tq by 
directly integrating Eqs. ^ and ([3]) at z = (solid line) and z = 2L£, (dashed line). 

of z/L]^ and t/tq based on Eq. f|T7j) . We see that the amplitude of the soliton undergoes 
only a slight decrease and its width undergoes a slight increase after propagating to a long 
distance. A simulation of the interaction between two ultraslow optical solitons is also 
carried out by initially inputting two identical solitons [see Fig. EJ^b)]. As time goes on, the 
two solitons collide, pass through, and depart from each other. They recover basically their 
initial waveforms after the collision. Finally, we have also made a numerical simulation by 
directly integrating Eqs. ([2]) and ([3]) to confirm the analytical prediction, with the result 
shown in Fig. [6](c). 
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B. A-type system 



The probe-field susceptibility for the A-type system is 



_ A4|ep-di2ppp31+P32 _ (1) (3)|p ,2 (22) 
Ap — Q — Xp ^ Xp I'^pl ! K^'^) 



where 



^(1) ^ ./Valep ■ ai2| / p2Pn_ P22_ 

(3) _ Ar,,|e,-di2|^ / + 5 + C -2A + i?*-C \ 

" eo^^ 1^ + j' ^23) 

with A = ipp'Sil/dl - l/rfi)/r, B = {ppf-^/di-ppf^/d* - i7]TA/2)/{2A), and C = [XiA - 
p\di\'^{d2B - d;B*) - \d2\\dlB - diB*)]/X2 with Xi = -iT{2\di\'^ - pl^ap), and X2 = 

ir(|(ii|2 + p|d2p). 

Shown in Fig. [5]^c) and Fig. [5]^d) are Re(x''^^) and Im(x^^)) as functions of the SGC 
parameter t]. The values of the system parameters are the same with those used for the 
V-type system. We see that both Re(x*''^'') and Im(x'^^)) decreases rapidly as r] grows, i.e., 
the SGC effect weakens the Kerr nonlinearity of the system. Based on this reason and on 
the large linear optical absorption shown in the last section, we conclude that an ultraslow 
optical soliton is not possible in the present A-type system. 



IV. THE CASE FOR OPEN SYSTEMS 



The results presented above are valid only for close systems. However, realistic physical 



systems, in particular molecules [35|], are usually open ones in which additional spontaneous 
decay pathways from upper levels to some other lower levels exist and hence spoil quantum 
coherence. In this section, we examine an open system related to V-type configuration and 
show that the additional decay pathways may lead to different initial population distribution 
of the states |2), and |3) and hence are detrimental to the SGC of the system. 

Consider an extended V-type system with its main part the same as Fig. [U^a), but the two 
closely spaced upper levels |2) and |3) may decay into the fourth state |4) with spontaneous 
decay rates and F3, respectively. In addition, a transient relaxation rate 7 is introduced 
to denote the leaving and entering of particles in light-particle interaction region; see Fig. [71 
Notice that in this simplified description the state |4) represents many other ground-state 
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FIG. 7: (color online) Model scheme of the open system related to V-type level configuration. Two 
closely spaced upper levels, |2) and |3), decay into the fourth state |4) with spontaneous decay 
rates T2 and respectively. \j) {j = 1, 2, 3, 4) are bare states, Ep is the probe laser field, A and 6 
are detunings. 7 is the transient relaxation rate of the particles leaving and entering light-particle 
interaction region. 



sublevels that provide various relaxation pathways from the two upper sates |2) and |3). 
The equations of motion for the density matrix of the system are 

Pu = -l{pu - Pll) + ^2p22 + TaPss - i^ppi2 + i%p2i - ip^ppis + ip^ppsi 

+?7A/f^(p23 + P32), (24a) 
P22 = -7(P22 - pll) - {T2 + r2)p22 + i^ppu - i%p2i - ri ^^f^ {p2z + P32), (24b) 



P33 = -7(P33 - P33) - (^3 + r3)p33 + ip^pPlZ - ip%Pzi - V ^^f^ (P23 + P32), (24c) 



P44 = -7(P44 - P44) + r2P22 + r3P33 

r2 + r'2 



2 (A + 5) - 



P21 

P31 

P32 = - ( «2A 



2(-A + 5) - 



2 

r3 + r'3 



P21 + i^pipil - P22) - ip^pP23 - V ^ P31, 



(24d) 
(24e) 



P31 + ipQpipu - P33) - i^pp32 - V — ^ — P21, (24f) 



12 + + r3 + , -o* ^ ■ o 

P32 - tilpPsi + W^lpPl-i 



-V n (P22 +P33), 



(24g) 



with p4j = (j = 1,2,3). Here p^^ (j = 1,2,3,4) is the population of thermal equilibrium 
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when the hght field and the spontaneous emission are absent. The equation of motion for 
the probe- field Rabi frequency Qp is still given by Eq. ([3]). Note that X]j=i Pjj 7^ 1 is broken 
for the present system. 

When the probe field is absent, the steady-state solution of the system reads 

, _ r(p22 + P33)+^r(p32+P23) , eq .o. x 

Pll — T Pii, [Zb&) 

P. - Mlf^±M^, (25b) 
P33 - MlP^M^^ (25c) 

(7 + r + r)(22A-r-r)r/r/2 

P32 = 7(P22 + P33), (25d) 

r'(p22 + P33) , eq /OK ^ 

P44 = h P44, (25e) 

7 

with other pji = 0, where Mi = 47(7 + T + T') + 7(r + V) [{T + T'f + 7(r + V) - rfV^ 12], 
M2 = 7(r + T')r]^Ty2, and M = -(T + T'){j + T + r)r]^T^ + (7 + T + T')^[(T + + 4A^]. 
For simplicity we have assumed r2 ~ = F and F2 ~ F3 = F'. One sees that there are 
non-zero initial population in the states |2), |3), and |4) due to non-zero Fg, F3, and pj^. 
It is easy to obtain the linear dispersion relation of the system 



c 



{UJ + fi'3)(Pll - P22) - ipif]T{pii - p33)/2 



D{uj) 

, P^{^ + ^2)(Pll - P33) - m^jpll - P22)/2 ] . . 

+ D{UJ) J' ^^^^ 

where D{uj) = (w + d2){oJ + 4) + 7]'^T'^/4: with c/g = A + 5 + i(F + F')/2 and c/g = -A + 
5 + i(V + F')/2. Illustrated in Fig. [Hl^a) and Fig. [Hl^b) are respectively the spectra of Re(i^) 
and \m{K) as functions of u for F' = O.OIF (the dashed line), F' = O.IF (the solid line), 



and F' = F 
be lO-^F 



^the dash-dotted line). When plotting the figure, the value of 7 is chosen to 
36| . In addition, rj is taken to unity in order to show clearly the infiuence of the 
decay rates and F3 with the maximum SGC From Eq. (126|) and Fig. [8] we obtain the 
following conclusions: (i) Due to the open character of the system (i.e. the existence of the 
additional spontaneous decay pathways from the states |2), |3) to the state |4)), the slope 
of Re(i^) curve with respect to u decreases when F' increases, which means that the group 
velocity of the probe field becomes larger when F' grows; see Fig[8](a). (ii) The absorption of 
the probe field can not be completely eliminated by means of the SGC effect, i.e. the system 
is always absorptive. In particular, the minimum of the absorption spectrum increases when 
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FIG. 8: (color online) Re(-fC) (panel (a) ) and lm{K) (panel (b) ) as functions oi uj/T for T' = O.OIF 
(dashed line), F' = O.IF (solid line), and F' = F (dash-dotted line). 



F' increases; see Fig|8[b). This tells us that the additional spontaneous decay pathways in 
the open system are detrimental to the SGC. Such conclusions can also be understood in 
another way. In order to suppress the absorption of the probe field, it is necessary to keep 
the most of atoms populated in the ground state |1). Nevertheless, from Eq. 1^5^ we see that 
the additional spontaneous decay pathways distribute the population of the ground state 
into other ground-state sublevels. In addition, when F' grows, the population in the state 
1 4) also grows, which results in the destruction of the quantum coherence of the system. 
Consequently, a quantum coherence with zero (linear) absorption based on the SGC takes 
place only in closed systems. However, as shown in Fig. El^b), although in an open system a 
high-quality quantum coherence can not be realized, a significant suppression of probe-field 
absorption can still occur, as shown by the recent experimental observation in a molecular 
system [35 1. 



V. DISCUSSION AND SUMMARY 



Notice that the SGC occurs in systems having near-degenerated levels with the same 
angular momentum quantum numbers J and mj and nonorthogonal dipole moments, which 
are rarely satisfied for realistic atomic systems 371]. However, such type of quantum inter- 
ference can be observed in many other systems such as semiconductor quantum wells and 



quantum dots |38l-l40| , autoionizing media [41[ , and anisotropic vacuum [42] . Our theoretical 



approach presented above can be easily generahzed to these systems with the SGC. 
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In summary, in this work we have investigated the hnear and nonhnear pulse propagations 
in hfetime-broadened three-state media with SGC. Three generic systems of V-, A-, and S- 
type level configurations have been considered and compared. We have shown that in the 
linear propagation regime the SGC in the V-type system can result in a significant change 
of dispersion and absorption and may be used to completely eliminate the absorption and 
largely reduce the group velocity of the probe field. However, the SGC has no effect on the 
dispersion and absorption of the A- and S-type systems. We have also shown that in the 
nonlinear propagation regime, the SGC displays different influences on Kerr nonlinearity for 
different systems. In particular, it can enhance the Kerr nonlinearity of the V-type system 
whereas weaken the Kerr nonlinearity of the A-type system. By exploiting the SGC effect, 
stable optical solitons with ultraslow propagating velocity and very low pump power can be 
produced in the V-type system by using only a single laser field in the system. 
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Appendix A: The expressions of A, B, and C in 



The expressions of A, B, and C in Xp and W are given as A — pXi + 2pX2 + X^, 
B^2Xi+X2+ pXs, and C = (2 + p)Xi + {2p + 1)X2 + pX^ + X*^A + B, where 



^3 = 



^1 = 



r(r2 + 4A2 - AC) 

2{Y2 - Y,)e + [r(y3 + n) + ^2A{Ys - Y,)]C - {T^ + 4A^)Y2 

r(r2 + 4A2 - AC) 

2{Ys - Y^)C + [(r + i2A)(yi + Y2)]C - r(r + 2iA)Ys 

r(r2 + 4A2 - AC) 



(Al) 
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with 



1^1 = 

l4 = 



0^3 - WC, ^ dl + ipC 



d2d3 + e 
ds - ipC 



dldl + 



+ 



d2d3 + 
p^d2 - ipC 



+ 



P^d* + ipC 
dldl + ' 
p^dl + ipC 



dl + ipC 

d2ds + e ' 



+ 



(A2) 



and C = ?7r/2. 
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